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Abstract 

The complex field zeros of the Random Energy Model are analytically determined. For T < Tc they 
are distributed in the whole complex plane with a density that decays very fast with the real component 
of H. For T > Tc a, region is found which is free of zeros and is associated to the paramagnetic phase 
of the system. This region is separated from the former by a line of zeros. The real axis is found to be 
enclosed in a cloud of zeros in the whole frozen phase of the system. 



PACS: 75.10.Nr (Spin Glasses), 05.20.-y (Stat.Mech.) 



1. Introduction 

The determination of the zeros of the partition function allows much insight into the behavior of 
a model. Since their relation with phase transitions were clarified by the works of Yang and Lee [1,2] 
and M. Fisher [3], many interesting properties of zeros distributions have been studied in homogeneous 
systems [4-9]. 

More limited, on the other hand, are the results concerning zeros distributions of disordered models. 
Taking into account that all the information about a given model is contained in the zeros of its partition 
function, it is reasonable to expect that many of the interesting properties [10] of disordered systems will 
be somehow reflected by them. 

An example of this can be found in the complex magnetic field zeros of spin glasses, which are 
dense near the real axis. This property has been suggested by Parisi [11] as being a consequence of the 
multiplicity of states and was recently illustrated [17] by a numerical calculation of the complex field 
zeros distribution in the Random Energy Model (REM) [12]. 

Here an analytical calculation of the complex field zeros of REM is reported. The technique is the 
same as recently employed by Derrida [13] to calculate the complex (3 zeros of the same model, which 
had been also numerically estimated [14,15]. 

In section 2 we describe the basics of the REM and give an account of the ideas that allow the 
calculation of its zeros. Section 3 is devoted to the REM in a complex field and the calculation of the 
dominant contribution to its partition function in each region of the complex h = (3H plane. In section 
4 these results are used to obtain the density of zeros of REM in the complex field plane, while section 5 
contains our conclusions. 

2. THE RANDOM ENERGY MODEL 

The REM is defined [12] to have 2^ independent random energy levels Ei, i=l,..2^ distributed 
according to P{E) — {nN)^^/^ exp{~E'^ /N). The partition function of a sample (a sample is a set of 
2^ numbers Ei), Z = X]i=i e~^^' can be written as Z = n{E) e"^^ where n{E) is the number of 
levels between E and E + A_E. This density of states n{E) if found to satisfy 



(n(£;))~exp7V(log2-(f)2) 

(1) 

{n{E)^) - {n{E)f ~ {n{E)) 

so two different energy regions can be identified: 

If \E\ < Ncc , witli €c = (log 2)^/^, then the number of states n{E) is much greater than one and its 
fluctuations are small, so for a typical sample we will have [13] 

n*yP(£) = {n{E)) + r]E{n{E)y/^ (2) 

where r/^ is a random number of order 1. 

On the other hand, if \E\ > Ncc , then the number of states is exponentially small meaning that a 
typical sample will have no levels in this region 

ntyP(E) = (3) 
The partition function for a typical sample will then be 

Z*yP = A + B (4) 

with 



|E|<JVec 

and 



B= Yl m{niE)y/^e-^^ (6) 

|E|<iVec 

The calculation of A may be done by the steepest descent method [16] for complex (3 — f3i + iji^- 

A - 2^y de exp -7V(e2 + /3e) (7) 

The integration contour is stretched to infinity through paths of steepest descent for the real part of 
the exponent of the integrand. When the integration limits lay on opposite sides of the saddle point in 
e = — /3/2 we can write 

(8) 



=/ +/ +/ 



and A has contributions both from the saddle point and from the integration limits for large N. This 
happens if |/3i| < Pc = 2ec = 2(log2)^/^ and we get 

A ~ exp{7V/3ee} + exp{7V(log2+(/3/2)2)} (9) 

whereas if > /3c the two steepest descent paths run from the integration limits down to infinity on 
the same side of the saddle point so it does not contribute and we have 

(10) 



= / + / 

-ec J —ec */±oo 



so 



A ~ exp{iV/3ec} (11) 

It has to be mentioned [13] that the lowest lying states in REM have order one fluctuations from sample 
to sample [12] so the contributions coming from the integration limits, both in (9) and (11), are known 
up to a fluctuating complex (if (i is complex) factor of order one. We will return to this point later on 
because it results crucial for our calculation. 

The contribution B of the fluctuations can be calculated for large N. Owing to the randomness of 
r]E,B is not the integral of an analytic function so the steepest descent method is no longer applicable. 
On the other hand, the r]E arc uncorrelated for different values of E so the term with the largest modulus 
will dominate the sum in (6). We can then estimate the modulus of B as 

\B\^ m^ \{n{E))y^e-^^ \ (12) 

—Nec<E<Nec 

It is easily seen that when \fii\ < /jc/2 the maximum is in between the limits so 

^2 

while if > /Jc/2 one of the limits dominates and 



B| ~exp{-(log2 + /3f)} (13) 



|B|~exp{7V|/3i|ee)} (14) 

We will now use these results to calculate Z for a typical sample of REM when a complex fleld H acts 
on it. We will take /? real from now on, as well as hi and non-negative [h = (3H = hi + ih^) which 
can be done without any loss of generality. 

3. REM IN A COMPLEX FIELD 

When a magnetic fleld H acts on the system, the energies of the states now have, in addition to 

the internal energies (coming from interactions among spins) a contribution from the interaction with 
the field, Ei = Ui — HM^, where the are the magnetizations of the states. The partition function 
Z = J2i e~^^' is then 



2" 



Z{l3,h)=y e-'5^'+''^' (15) 



i=l 



The internal energies Ui have a Gaussian distribution which is independent of the field, so taking into 
account that among the 2^ states of the system, = Qm = ( (w+m) ) have magnetization M = Nm, we 

2 

can write 

Z{p,h) = Yl e^'"''^-(/3) (16) 

-l<m<l 

where 

ZM=f^ e-^""' (17) 

is the partition function of a sample of REM with gm instead of 2^ states. If we write gm = ((Zm)^, then 
all the results of last section concerning Z are valid for Zm with the sole condition of replacing log 2 by 

loggm = log2 - i [(1 + to) log(l + to) + (1 - m) log(l - to)] 



(for large N). 

These restricted partition functions can be written as 



2m{P)= Yl ri{u,m) e-"""^ (18) 

— 00<li<00 

where n{u, m) are the numbers of states with internal energy Nu and magnetization Nm. They can be 
seen to satisfy 

(n{u,m)) ^ cxp Af(loggf„ - u^) 

(19) 

{n{u,ni)'^) — {n{u,m)y ~ {n{u,m)) 
Then for a typical sample we will have 

{n^^P{u,rn) ^ {n{u,rn)) + ri{u,m){n{u,m)y^^ ,i{ \u\ < Uc{m) ^^^^ 
n*yP(u,m) ,if |u| > Uc{m) 

where the ri{u, m) are random numbers describing the typical fluctuations, and the energy threshold 
Uc{m) is now m-dependent. 

Ue(m) = (logg„)V2 (21) 
The partition function of a typical sample under field will then be 

2*yP(/3, h) = Y, e'^'^^Z^iP) (22) 

m 

We now use the results of last section and write 



Z^yP{(3)=Arn{P)+BM (23) 



with 



and 



An = Pm exp {NPuc{m)} , if /3 > 2uc{m) 

Am = Pm exp {NI3uc{m)} + exp {Ar(log g„ + {13 /2f)} , if (3 < 2uc{m) 

Bm = Vm exp {NPuc{m)} ,if (3 > Uc{m) 

Bm = Vni exp {^i^oggm + f3^)} ,if (3<Uc(m) 



(24) 



(25) 



Here "qrn is a random number of order one which describes the fluctuations in the density of states as 
already mentioned. The origin of pm is somewhat different. They arc due to fluctuations in the values of 
the lowest lying states: The ground state energy (internal energy in this case) Uq has order one sample 
to sample fluctuations [12,13] so U^'^im) ~ N{uc(m) + Q/N) with C, an order one random number so 
ef^^° ~ pe^^"'=^'"\ where p is a positive fluctuating number. 

We see in this way that both in A and B the ground state contributions are non-analytic (due to 
disorder), and therefore non-integrable by the steepest descent method. The only analytic contribution 
is the saddle-point term in (24). 

Let us now turn to the evaluation of the sum in (22). We define (3c{m) as 



/3c (m) = 2uc{m) 



(26) 



and its inverse rndP) such that satisfies 



(3/2 



log 2 — - [(1 + TOc) log(l + mc) + (1 — mc) log(l — mc) 



1/2 



(27) 



The relation hdP) — atanh(mc(/3)) defines the critical field he which is the limit between the frozen and 
paramagnetic phases of REM with field [12]. This critical field satisfies 



/3/2 = [log 2 + log(cosh(/ic)) - /ictanh(/ic)]^/^ 



(28) 



Using (26) wc can rewrite the condition /3 > 2uc{m) as jmj > niciP) and the condition (3 > udm) as 
\m\ > mc(2/3). The sum in (22) can now be written 



z^yp = A + B 



(29) 



with 



Nmh 



and 



Nmh 



We then have 



(30) 



(31) 



A = Ai + A2= 



Pme 



„iV[(/3/2)2+logg„.+m/t] 



-1<TO<1 



|m|<mc(/3) 



(32) 



B = Bi + B2 



|jn|>mc(2/3) 



rj^e^l0'+^-s,rn+2mh] (33) 



The only analytic contribution is A2, which will be estimated by means of the steepest descent method. 
Wc will ignore Bi because it is already included in Ai . Let us first discuss the evaluation of the fluctuating 
contributions Ai and B2. The non-analyticity of the integrand makes it impossible to use the steepest 
descent method. On the other hand, as already discussed in section II, the term with the greatest modulus 
will dominate the sums so we have 



max exp{N[P{\ogqm)'/^ + mhi]} (34) 

— l<m<l 



and 



B2I ~ niax exp {^- [13'^ {log qm) + 2mhi]} (35) 

|TO|<mc(2/3) 2 



It is worth noticing that 



\B2\=0 if/3>/3c/2 



(36) 



On the other hand, if /3 < /3c/2,then mc{2(3) > and there exist two ranges of hi in which B2 gives 
different contributions: 

The maximum of |i?2| is at m = tanh(2/ii) when \ fh\ < mc{2(3) and at ±mc(2/3) when |m| > mc(2/3). 
This can be written, using (27) and (28) as 

rexp{f [/?2 + log2 + log(cosh(2/ii))]} ,if \hx\ < \K{2(i) 
\B2\^{ (37) 
[exp{iV[/32 + |/ii|m^(2/3)]} ,if I > 1/1,(2/3) 

Now we discuss the evaluation of Ai . The maximization condition is satisfied at m = tanh(/i) where h is 
a function of hi and /? that satisfies 



^h_ 
2h'i 



log 2 + log cosh(/i) — h tanh(/i) 



1/2 



(38) 



If we put h = hc{P) with /3 = a/3 , now the unknown is a(/3, ft.i). Replacing this in (38) we find h ■ 
so the following set of equations: 



ahi 



1 /3/2 = [log2 + logcosh(/i) - /itanh(/i)]V2 ^^^^ 
1 /3/i = /3/i 

is equivalent to (38). 

It is easy to see that when hi = hc{(3) wc obtain a = 1 while a < 1 \i hi > hc{f3) and a > 1 if 
hi < hc{P), so h will always be between hi and hc{P). 
The term Ai will then have the following form 

|Ai| -cxp{A^[^ + /),itanh(/i)]} (40) 

The calculation of A2 may be done by the steepest descent method (sec the appendix for details) for 
large A'^. It follows that again two different expressions arc found, now depending on the values of h and 
(3. A smooth arc is found in the complex h plane which divides these two possibilities. This arc touches 
the real axis at hc{P) and the imaginary axis at h = in/2. In the inner part of this arc we find 



A2 ~ expiV{^ + hmdP)} + e^pN{^^^ + log(cosh(/i))} (41) 



and in the outer side 



A2~exp7V{^+/iTOe(/3)} (42) 

It is clear that Z will in general be the sum of several terms of the form e^*^ with ip a complex function, 
so for large N just that term with the greatest modulus will be relevant. We then have to determine, in 
each region of the complex h plane, the term with the greatest real part of if. This will give rise to the 
appearance of a phase diagram for the model in complex h. 

In what follows we will be only interested in the determination of the modulus of Z, which, as we 
shall see, is enough to calculate the density of zeros of the model. If we write \ Z \= expNcj), the different 
contributions we have to compare are 



= ^!-±^ + log I cosh(/i) I for /3 < /?c and I /ii I < h}}"^ (/i2 , /?) 
<^2 = Y + |/ii|mc(/3) for /3</3e and V/ii 

03 = ^ + /ii tanh(/i) V/3 and V/ii 

04 = ? + ^ + ^log(cosh(2/ii)) for and \hi\ < lhc{2(3) 



05 = + |/ii|tanh(/ie(2/?)) for /3 < ^ and \hi\ > ^/ic(2/3) 

Three different ranges of temperature can be identified; 

a) ff /3 > 00 then only 03 exists so = 03 , V/i. 

b) If /3c > /3 > /3c/2 we get contributions from 0i, 02 and 03. In this case two different regions are found 
in the complex h plane (figure 1, upper portion), according to what term is the most important. These 
regions are separated by an arc which goes from (/?,c(/3), 0) on the real axis to (0, /i2(/3)) on the imaginary 
axis. In the inside of this arc 0i dominates while in the outside it is 03 which is the most important. 02 
is never relevant. 

c) If /3 < (3c/2 then also (p^ and 05 contribute, and they have to be compared to 0i and 03 in their 
respective zones. It is found that 05 is never relevant, i.e. (0i — 05) and (03 — 05) are never negative. 
On the other hand 04 is greater than 03 (for hi < ^hc(2/3) ) and it is also greater than 0i in the upper 
portion of the zone where 0i > 03, so the line where 04 equals 0i is lower than that where 03 equals 0i. 

The situation in case c) is that depicted in the lower portion of figure 1. The arc separating 0i and 
04 now touches the imaginary axis at a point (0, h^i^)) which results from equating 04 to 0i for hi = 
and satisfies 



^u=v'^2(/?)) =e"(^^^ (43) 



so we can see that h^iP = 0) = 7r/4. 



4. DENSITY OF COMPLEX FIELD ZEROS OF REM 

Having obtained = '"^'^ in the various phases of the model, we can calculate the density of zeros p{h) 
of Z by means of the formula [13] 

1 1 / /)2 a2 \ 

p(/,)= V20=- + — 2 0(/ii,/i2) (44) 



27r ' 27r \dhl dhl 

which is a consequence of the electrostatic analogy [2] that identifies p{h) with a charge density and 
with the resulting electrostatic potential. 

The density of zeros associated to 0i is null (because the real part of an analytic function has a 
vanishing Laplacian). The densities of zeros corresponding to 03 and 04 are respectively 



P3 = ^VV3 = a(l - tanh^(a/ii)) ^\ (45) 

27r + 2h{{l - tanh (a/ii)) 

Pi = ;^VV4 = 2(1 - tanh=^(2/ii)) (46) 



where a = h/h\ as already defined. 



In addition to these distributed zeros, there exist other zeros located right on the Unes dividing the 
different zones, and which density is proportional to the discontinuity in the normal component of the 
"electric field" across that line. On a line separating two generic phases a and b we then have 



Pf = ^\^i<f>a-M\ (47) 



The density of zeros on line 1-3 is found to be 



13 1 r/ 1 /?\ sinh(2/ii) 2 / sinh(2/ii) 2i /.q\ 

Pl = — |(tanh(/i) - cosh(2/zi)+cos(2/i2)^ + ^cosh(2/zi) + cos(2/i2) ^ ^ ^^^> 

On line 1-4 the expression is the same with the only difference that h{h = ahi) is replaced by 2hi. Their 
zero-densities are equal where lines 1-3 and 1-4 meet. 

The line 3-4, which is the vertical at hi = ^ft,c(2/3) is not a line of zeros, i.e. the electric field is found 
to be continuous there. On the other hand the density of distributed zeros jumps as one crosses line 3-4 
by a factor {1 + 2{hi/l3f{l - tanh^(2/ii))}. 

The density of zeros on the line 1-3 is found to be null at the point where that line touches the real 
h axis, in agreement with the second order character of the field-driven transition of the model. 

5. CONCLUSIONS 

We have estimated the behavior of the partition function of REM in complex field. Three different phases 
appear: 

For T <Tc there is only one phase which we called phase 3 and is the continuation to complex field 
of the frozen phase of the model. For T,, < T the paramagnetic phase is also present, which we call phase 
1 (figure 1, upper portion). Finally for T > 2Tc what we called phase 4 also appears at complex values 
of h (figure 1, lower portion). There is no physical (i.e. for reaZ values of h) counterpart for this phase. 

We have also calculated the densities of zeros in the different phases. A particular structure has 
been revealed: The zeros are dense near the real axis in the frozen phase. This characteristic of the zeros 
distributions is intimately related to the properties of spin glasses and has been discussed by Parisi [11], 
although no explicit calculation had been done before. A numerical estimate of the complex field zeros 
of REM has been done [17] and shows good coincidence [18] with the results here reported. 

The existence of distributed zeros is in this model related to disorder in a very clear way: The 
fluctuating factors pm and rjm are responsible for the appearance, in Z, of nonanalyticities of the sort 
which are necessary to obtain a non-zero Laplacian and consequently distributed zeros. 
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APPENDIX 

Here we describe the evaluation of A2 by means of the steepest descent method [16]. Starting from 

A2= dm e^*(™) {Al) 

with 

*(m) =mh + log2 - ^ [(1 + m) log(l + m) + {l-m) log(l - m)] + ^0 (^2) 
we see that A2 is non-zero only for (3 < Pc, otherwise mdP) = and this term does not contribute. 



Following the usual procedure for evaluating this kind of integrals in the complex plane, we stretch the 
integration contour to infinity following paths of steepest descent for the real part of ^ . The imaginary 
part of the integrand will in this way be constant along the whole path, and may be taken out from the 
integral. The level curves of ^ are topologically similar to those of z^. Two steepest descent lines cross 
each other at the saddle point which is located at mo = tanh(/i): one of them goes down to the valleys to 
the left and right of the saddle-point and will be called AA' while the other climbs towards the mountain 
tops and is the line LL'. The way in which the original integration path has to be deformed depends on 
the relative positions of ±mc and LL'. 
Depending on (3 and h, two situations are possible: 

a) When the integration limits ±mc(/3) lay on opposite sides of LL', there are contributions from 
three paths: the one that runs through the saddle-point (path AA'), and those ending at ±mc(/3). 

b) When both integration limits are on the same side of LL', the path that joins the valleys through 
the saddle-point does not contribute to the integral. 

The limit condition for this two cases is that mc(/3) lays exactly on LL', which in turn is equivalent to 
the condition that Imag(\l/(mc)) =Imag(\E'(mo)). This last identity results from the fact that LL', being 
a steepest descent for the real part of is a level curve for the imaginary part of ^f. This limit condition 
is then written 

tanh(/^0 = ^ ,^^^/^lim) ^^3^ 

tg(n2) 

and for /3 < /3c it defines a smooth arc in the complex h plane (with hi and /12 positive) that joins the 
points {hc{(3), 0) and (0, 7r/2). The shape and size of this arc depend on (3. In its inner part situation a) 
is valid so for large N 

A2 ~ exp A^{^ + hnieiP)} + exp7V{^-^^ + log(cosh(/i))} (A4) 
while in the outer side of this arc b) holds and 

A2~exp7V{^+/ime(/3)} (^5) 
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FIGURE CAPTIONS 

FIGURE 1. Shown are the phases of the model in the complex field plane for two different values of the 
temperature. For < T < 2Tc(upper figure) two distinct zones appear separated by an arc of zeros. 
Above this arc (phase 3) densely distributed zeros are found while the region below it (phase 1) is free 
of zeros. For T > 2Tc the outer part of the arc appears divided in two zones by a vertical line, which is 
not a line of zeros. The densities of zeros are different to the left (phase 4) and right (phase 3) of this 
vertical line. 
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